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A SPECIAL  DISTRIBUTIONAL  RESULT  FOR  BILINEAR  FORMS 


Janes  Beckett  III,  William  R.  Schucany,  and  James  D.  Broffitt* 


ABSTRACT 


Necessary  and  sufficient  conditions  are  given  such  that  a quadratic  form 
has  moment  generating  function  E(exp(ty'  B U)]  - (l-t2)”r/^  for  |t|  < 1 with 
U - \ (y,  f)  *nd  £ positive  definite.  An  important  corollary  gives  conditions 
under  which  the  bilinear  form  X'AY  involving  two  different  multivariate  normal 
random  vectors  (of  not  necessarily  the  same  dimensions)  has  the  same  distribution 
as  the  sum  of  independent  random  variables,  each  having  the  LaPlace  (double 
exponential)  distribution. 


KEY  WORDS:  LaPlace,  double  exponential,  correlation.  Normal,  products 
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1.  MAIN  RESULT 


Let  Z be  a random  variable  which  has  the  same  distribution  as  one  half  the 
difference  of  two  independent  x2(v)  random  variables.  Then  the  moment  generating 
function  of  Z is  M^ft)  • (l-t*-)  V^2  for  |t|  < 1.  We  will  use  the  notation 
Z * L(v)  to  denote  this  fact.  There  is  an  interesting  relation  between  the 
distribution  of  Z and  the  LaPlace  or  double  exponential  distribution.  If  W 

has  a LaPlace  distribution,  i.e.,  W has  probability  density  function 

1 - Iw  1 2 -1  i i 

f (w)  - -z  e ' ' then  the  moment  generating  function  of  W is  (1-t  ) for  |t|  < 1. 

2 D v/2 
Thus  if  v is  even  and  W , . . . »wv /2  are  ****  LaPlace  variates,  then  Z *-  E W . 

VD  m , i“l 

If  v is  odd  (v-2m+l),  then  Z - l Wt  + (Sj-S2)  where  S,,  S2 

are  mutually  independent,  the  W's  are  LaPlace  variates,  and  and  S2  are 

chi-square  variates  each  with  one  degree  of  freedom.  The  sum  of  LaPlace 

variates  characterization  for  the  distribution  of  Z is  convenient  for  calculating 

probabilities.  For  example,  if  c > 0,  then 

Fjfc)  - P(WX  < C)  - 1 - I e"c  ; 

F2(c)  - P(WX  + W2  < c)  - F1(c)  - f e_C  ; 

F3(c)  - P(W j + W2  + W3  < c)  - F2(c)  - e“C  . 

These  results  were  developed  from  the  general  result  for  the  p.d.f.  of  the 

average  of  LaPlace  variates  given  in  Johnson  and  Kotz  (1970,  p.  24).  In  fact 
m 

the  p.d.f.  of  I W is  given  by 


i-1 


- w 


f(v)  - 


V 2J(2m-l-2)l  |w|j 

* 1 / ■ « \ • • 


22m"1(m-l)!  j-0 


The  theorem  and  corollary  presented  below  establish  necessary  and  sufficient 

conditions  for  the  bilinear  form  in  normal  variates,  X'AY  to  be  distributed  as 

L(v).  This  is  followed  by  some  examples  one  of  which  provides  a new  exact  test 

for  the  correlation  coefficient  when  Oj  and  o2  are  known  by  giving  the  (null) 

n 

unconditional  distribution  of  r (X .-X) (Y .-Y) . Also  a known  result  [Nyguist 

i-1  1 


et.  al.  (1954),  Nicholson  (1958),  Mantel  (1973)]  is  demonstrated  as  being  a 
special  case  of  the  corollary. 


Theorem;  Assume  U * N^  (jj,  f)  with  I positive  definite;  then  E(eC'  ^)  • (1-t2) 
for  -1  < t < 1 if  and  only  if  rank  (B)  • r.  Bp  - 0,  tr  (Bf)  ■ 0,  and  4BZBI 
is  idempotent. 


Proof:  The  quadratic  form  can  be  rewritten  in  the  following  manner: 

u* bu  - u’z"1/2(z1/2bz1/2)  z“1/2  y - u’  z“1/2  (9  px  9')  z“1/2  u 

- 2’  Dx  Z 

k 2 

“ 1 Ai  Zi 

i-1  1 1 

1/2  1/2 

where  Q Dx  Q'  ■ I B I , § Is  orthogonal,  X's  are  eigenvalues  of 

Z1/2  B Z1/2,  and  Z - Q*  Z~1/2  U * N (Q'z“1/2  p - 0,  I).  Then  the  m.g.f. 

k 

of  U* BU  can  be  expressed  as  M.,,B1I(t)  - n M_2  (tX.) 

^ K i-1  Zi  1 

2 

k 0*tX  Mi 

(1.1)  • n (exp  lyijdi1*  <1“2txi)"  ' ) if  tXi  < 1/2  for  1-1,..., k. 

Note:  (i)  Bp  - 0 «■— > Dx©  - 0 <—■>  0^  - 0,  i-l,...,k. 

k 

(ii)  tr (Bf ) - Z X.. 

i-1  1 

(iii)  4BZBZ  idempotent  <— > 4 Z1/2  BZBZ1/2  idempotent. 

(iv)  4f 1/2  »IK1/2  - QQU2  Q'  . 

(<— ) Now  assume  rank(B)  - r,  Bp  - 0,  tr(BZ)  ■ 0,  and  4BZBZ  is  idempotent. 
From  (1.1)  and  (1)  one  can  obtain 


VBU  <c>  - n (l-2tX.)‘1/2  , tX  < 1/2  i-1,...,  k. 

i-1 

Also  rank  (4Z1/2BZBZ1/2)  - rank(B)  - r . 


1/2  l/’ 

Thus  from  (iv)  and  the  idempotency  of  4E  BIBE  ” (by  (HO)  we  have 
2 

4X  "1  with  multiplicity  r 

- 0 with  multiplicity  k-r. 
k 

This  together  with  Z X^  - 0 (see  (ii))  implies 

i-1 

X - 1/2  with  multiplicity  r/2 
—1/2  with  multiplicity  r/2 
■ 0 with  multiplicity  k-r 

/.  My'BU(t>  - (l-t)“r/A  (l+t)-r/4  , | t < 1/2  and  - | t < 1/2 

- (l-t2)'r/4  , -1  < t < 1 . 


(“■>  ) Now  assume 

My.g^t)  - <l-t2)"rM  , -1  < t < 1 


Then  from  (1.1) 


k ^ -1/2  2 -r/4 

exp  { Z ) * J1  (l-2tX.)  1/2  = (1-t2)  X,U 

1-1  i-i 


-1  < t < 1 


k 0 tX 

(1.2)  /.  exp  {-4  Z ) 

i-1  wtAi 


n ( 1— 2t X . ) 2 = (l-t2)r  , — < t < 

i-1 


The  r.h.s.  of  (1.2)  is  a polynomial  of  degree  2r  with  2r  real  roots. 


(r  roots  - +1,  r roots  - -1),  thus  so  is  the  l.h.s.  But  roots  of  the 

^ 2 

l.h.s.  must  be  roots  of  n (l-2tX  ) . Thus 

i-1  1 


X - 1/2  with  multiplicity  r/2 
(1.3)  —1/2  with  multiplicity  r/2 

- 0 with  multiplicity  k-r 


and 


k o:t» 

exp  { -4  E 11 


i-1  1_2tXi 


} . 1 by  (1.2),  or  equivalently 


(1.4) 


k 0‘tA 


. 0 . 


From  (1.3)  and  (1.4)  it  la  easily  deduced  that  0^  • 0 If  j*  0, 
so  that  • 0,  1*1,...,  k.  Then  by  (1)  we  have  B^j  ■ 0 . 

From  (1.3)  and  by  (11)  we  have  tr(BE)  • 0.  Also  by  (1.3)  rank(E1^‘Bf^“)-r 
since  there  are  r nonzero  roots,  thus  rank(fc)»r. 


Finally,  by  (iv)  and  (1.3) 


4E1/2BEBE1/2  - Q 


(Hr) 


where  g " I9i 
kxr 


§•>) 


But  gjQ*  is  idempotent  since  Q ' - } . 

Thus  by  (ill)  4BEBE  is  idempotent.  Q.E.D.  So  under  these  necessary 
and  sufficient  conditions,  U'BU  * L(r/2). 


Corollary:  Given 


with  E positive  definite  and  the  (p+q)  x (p+q)  matrix  C defined  by 
C » f pxq  | , then 


/?  * \ 
“ f pxq  I 

I V 0 / 


qxp 

E [e£  * "]  - (1-t2)  v^2  for  -1  < t < 1 if  and  only  if  rank  (A)  • v, 
A Hj  ■ 0,  A'jfj  ■ 0,  tr  (AEj^)  * 0,  and  CJQ?  is  idempotent. 


Proof:  Let  B • (1/2)C,  then  $'*Y  - y'BU.  Thus  E(et'&y)  - (l-t2)"2v/4,  -1  < t < 1 

<— > (by  theorem) 

rank(B)  • 2v,  By  - 0,  tr(B£)  ■ 0 and  4$£B£  idempotent.  But  it  is  easily 
seen  that  these  conditions  are  equivalent  to  rank (A)  ■ v,  Ay2  ■ Q, 

A'y^  ■ 0,  trtAIjj)  • 0,  and  CICE  idempotent.  Q.E.D. 

2.  EXAMPLE  APPLICATIONS  AND  DISCUSSION 


Example  1.  If  we  consider  2p  independent  standard  normal  variates  X.,  X.. , ....  X 

1 * P 

and  Yj,  Y2,  ....  Yp  (hence  * ?21  * -p*  w'iere  tp  denotes  a pxp  zero  matrix, 
A • Ejj  • I22  • Ip,  and  ■ W2  * >> • then  " E XiY1  - L(p). 


Example  2.  If  X and  Y are  p-dimensional  and 

OH  ;))• 


where  Jp  denotes  a pxp  matrix  of  ones. 


then  X’(Ip  - ~ Jp)  Y - ^ (Xt  - X)(Y1  - Y)  « L(p-l).  The  matrix 

A - I - — J 

~p  p *p  is  a familiar  idempotent  matrix  of  rank  p-1  which  is  clearly 
orthogonal  to  the  vector  of  ones;  all  conditions  of  the  corollary  are  met 
by  Inspection. 


Example  3.  If  X^,  Xj,  X^,  X^  are  lid  standard  normal  random  variables,  then 

XjXj  i X3XA  has  a Laplace  distribution.  Letting  - (Xj.Xj),  *'-(X3,X4), 
the  result  follows  directly  from  Example  1.  The  problem  of  the  distribution 

|x i hi 

of  the  random  normal  determinant  D *|  X^  X2  | , was  solved  by  Nyquist 
et  al.  (1954)  and  generalized  by  Nicholson  (1958),  although  the  form  of  the 
p.d.f.  in  Nicholson's  result  is  off  by  a scale  factor.  [See  also  Johnson 
and  Kotz  (1970),  p.  30.]  Mantel  (1973)  used  conditional  expectations  and 


characteristic  functions  to  show  the  distribution  of  D - X.X.  + X0X . as 

1 > i U 

Laplace.  Clearly  in  this  instance  X^,  + X^X^  and  X^2  - X^X^  have  the 
sane  distribution. 


Example  4.  In  considering  the  asymptotic  distribution  of  an  inner  product  of 

independent  vectors  of  Friedman  rank  totals  Hollander  and  Sethuraman  (1978) 
recently  presented  a result  for  which  one  of  their  special  cases  bears  an 
interesting  relation  to  present  result.  In  the  usual  setting  the 
rankings  (1,...,  k)  in  each  of  m rows  are  summed  to  obtain  rank  totals 
Sj(j»l,...,  k)  and  from  the  multivariate  central  limit  theorem  one  concludes 
that  the  large  sample  distribution  of  S is  multivariate  normal  (singular 
since  ESj  ■ mk(k+l)/2).  If  we  have  a second  independent  table  of  n rankings 

with  totals  denoted  by  Tj  and  consider  the  case  of  random  rankings,  then 

2 1/2 

the  standardized  totals  X^  - [S^  - m(k+l)/2 J/fmk(k  -1)/12]  and 

Y " lTj  " n<k+D/2l/tnk(k2-1)/12]1/2,  j-1 k-1  are  as  m,  n - ~ (2.1) 

asymptotically  jointly  normal  with  means  ■ y2  ■ 0 and  (k-1)  x (k-1) 

* **  X 

covariance  matrices  ■ Z22  - (k-1)  (kl^  ” -k-1^' 

Hollander  and  Sethuraman  (1978)  find  the  asymptotic  distribution  of  a 

standardized  inner  product,  namely 

H - 12  {S’T  - mnk(k+l)2/4)/{k(k+l)(mn)1/2) 

to  be  that  of  ZV,  where  Z and  V are  independent,  Z is  standard  normal  and 
2 2 

V •v  x with  k-1  degrees  of  freedom.  It  may  also  be  seen  that 


-1 


H - X’  A Y 

with  X and  Y as  given  in  (2.1)  and  A ■ (k-l)k  x tlk_j  + J ^ ] . 

It  is  easily  verified  that  rank(A)  ■ k-1,  Ay2  - A * y ^ ■ 0, 

£ 1 2 ■ £21  " 0 *“>  tr(AE2^)  " °*  and  -'ll  " --22  ” ^k-1  ----  idemPotent 

and  therefore  the  limiting  distribution  of  H is  L(k-l).  This  distribution 

may  be  more  desirable  in  a setting  which  requires  the  calculation  of  quantiles 


mmmm 


mmmmmmm 


or  tail  areas.  Nevertheless,  we  note  the  interesting  distributional  result 
that  ZV  is  distributed  as  the  sum  of  (k-l)/2  independent  Laplace  variates 

C2  2 \ 1/2 

^3  + *4  / 

are  identically  distributed. 

Examples  1 and  2 essentially  provide  an  exact  test  for  the  correlation 
coefficient  p • 0 with  and  Uj  known  or  unknown  for  an  even  or  odd  number 
of  bivariate  normal  observations  with  known  variances. 
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